Electrodynamics of spin currents in superconductors 
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In recent work we formulated a new set of electrodynamic equations for superconductors as an 
alternative to the conventional London equations, compatible with the prediction of the theory of 
hole superconductivity that superconductors expel negative charge from the interior towards the 
surface. Charge expulsion results in a macroscopically inhomogeneous charge distribution and an 
electric field in the interior, and because of this a spin current is expected to exist. Furthermore, 
we have recently shown that a dynamical explanation of the Meissner effect in superconductors 
leads to the prediction that a spontaneous spin current exists near the surface of superconductors 
(spin Meissner effect). In this paper we extend the electrodynamic equations proposed earlier for 
the charge density and charge current to describe also the space and time dependence of the spin 
density and spin current. This allows us to determine the magnitude of the expelled negative charge 
and interior electric field as well as of the spin current in terms of other measurable properties of 
superconductors. We also provide a 'geometric' interpretation of the difference between type I and 
type II superconductors, discuss the relationship between our model and Slater's seminal work on 
superconductivity, and discuss the magnitude of the expected novel effects for elemental and other 
superconductors. 

PACS numbers: 



I. INTRODUCTION 

The theory of hole superconductivity predicts that su- 
perconductors expel negative charge from the interior to- 
wards the surface 0, 0, HI, resulting in the existence of 
an electric field in the interior of superconductors. To 
account for this fact we have proposed a new set of elec- 
trodynamic equations for superconductors 0, H| in place 
of the conventional London equations Q . The new equa- 
tions predict a charge density p(f) in the interior of su- 
perconductors that satisfies the differential equation (in 
the static case) 



V 2 p(r) = - T (p(r)- Po ) 



(1) 



with po & positive constant. This equation (together with 
a similar equation for the electric potential) predicts that 
the charge density in the deep interior of a superconduct- 
ing body is po, and that an excess of negative charge den- 
sity p- exists within a London penetration depth (Xl) of 
the surface. It also gives rise to an electric field in the in- 
terior of superconductors that points towards the surface. 
For a cylindrical geometry and cylinder radius R » Xl, 
charge conservation implies that 



Po 



2Xi 



P- 



and for a sphere of radius R >> Xl 

3A L 



Po 



R 



-P- 



(2a) 



(2b) 



In both cases, the electric field increases linearly away 
from the center, attains its maximum value 



-AttXlP- 



(3) 



at a distance Xl from the surface, and drops to zero at 
the surface. 

Energetic arguments Q indicate that p_ and E m 
should be independent of the size of the sample and that 
E m should be related to the square root of the supercon- 
ducting condensation energy per unit volume. However, 
the precise value of could not be determined from the 
treatment of refs. [H, Gj • The determination of the value 
of E m is one of the central results of this paper. The 
interior charge density po depends on the sample dimen- 
sions (Eq. (2)) and becomes smaller as the sample size 
increases. 

We have also pointed out in earlier work that an elec- 
tric field in the interior of a superconductor should give 
rise to a macroscopic spin current 0, @| , and showed that 
in fact the microscopic Hamiltonian proposed by the the- 
ory of hole superconductivity^ favors the existence of 
such a current [6| . In very recent work we have put forth 
a more detailed picture of how the charge expulsion pro- 
posed to exist in superconductors occurs and how the 
spin current is generated 8]: namely, that electrons 'ex- 
pand' their wavefunction|9] from being confined to a lat- 
tice spacing (as corresponds to a nearly filled band) to 
a much larger (mesoscopic) extent. As electrons move 
radially outward, the spin-orbit interaction gives rise to 
azimuthal velocities resulting in the two members of the 
Cooper pair circulating in opposite directions in an orbit 
of radius 2Xl, which encloses precisely one flux quan- 
tum of spin- orbit flux. This picture also provides for a 
dynamical explanation of the Meissner effect 0]. It pre- 
dicts that in the presence of an external magnetic field, 
a charge and a spin current will circulate near the sur- 
face of a superconductor. In the absence of an external 
magnetic field, a pure spin current will circulate, with a 
universal expression for the magnitude of the spin current 



2 



speed at the surface v® = 7l/(4m e Ai,).[l] 

For the charge current and charge density we derived 
their space and time dependence in Ref. [4| (the charge 
current behavior predicted by our equations is the same 
as predicted by the conventional London equations). In 
this paper we derive the equations governing the space 
and time dependence of the spin current, and show how 
the spontaneous spin current is related to the expelled 
negative charge. Furthermore we provide a new interpre- 
tation of type I and type II regimes in superconductors, 
and we explain how superconductors manage to conserve 
angular momentum. We also calculate the magnitude of 
the effects predicted for various materials. Finally, we 
discuss the relation of our theory with earlier pre-BCS 
work. 

We should point out at the outset that the aim of this 
paper is limited to providing a macroscopic description of 
the electrodynamics of charge and spin currents in super- 
conductors, and is not to provide a full microscopic de- 
scription of the superconducting state. Just like conven- 
tional BCS theory is consistent with conventional London 
electrodynamics, a full quantum-mechanical description 
of the problem should be compatible with the electrody- 
namics discussed here and will be the subject of future 
work. 



II. SPIN CURRENT CONSTITUTIVE 
RELATION 



We have shown in ref.[8j that the Meissner effect can 
be understood 'dynamically' by assuming that electrons 
move radially outward a distance 2Az, in the presence of 
an unscreened magnetic field B that gives rise to a vector 
potential 



A = 



B x r 



(9) 



and acquire through the action of the Lorentz force an 
azimuthal velocity given by Eq. (8). In a cylindrical 
geometry, the magnetic field at the surface equals the 
applied magnetic field B, and the vector potential at the 
surface is given by 



A(R) = X L B x h 



(10) 



(n = f denotes the direction normal to the surface) which 
coincides with Eq. (9) for r = 2Al- Eq. (10) results from 
solving London's equation for an infinitely long cylinder 
of radius R >> Ar. [lOfl, or from simply assuming that the 
magnetic field penetrates a distance Al and using that 



A - dl = B-da 



Similarly, we showed in ref.[8| that an 'effective' 
unscreened magnetic field B a acts on the supernuid elec- 
trons arising from the spin-orbit interaction 



B n 



2Tm s fl 



(11) 



For the charge current 

J = en s v s 



(4) 



with n s the superfluid density and v s the superfluid 
charge velocity, the governing equation is 



with 



eh 
2m P c 



(12) 



the electron magnetic moment. Electrons acquire a spin 
current velocity[8j 



J(r) 



4^A| 



(5) 



with A the magnetic vector potential, related to the mag- 
netic field B by 

VxA = B. (6) 
The London penetration depth is given by 
1 47m ,e 2 



\ 2 L m e c 2 



and Eq. (5) is equivalent to 



v s (r) = A(r) 

m e c 



(7) 



(8) 



which can be understood as resulting from the 'rigidity' 
of the wave function that 'forces' p in the relation p = 
Tie^s + {e/c)A to stay zero at all times and locations in 
a simply connected superconductor. 
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4m e Ai 



a x r 



-v'a x r 



(13) 



through the action of the 'Lorentz' force from the field 
Eq. (11) in moving radially outward a distance 2Al. 
The spin-orbit vector potential corresponding to 'mag- 
netic field' Eq. (11) at radial distance 2Al is, from Eqs. 
(11) and (9) 



A n 



2im s \Lfl x f. 



(14) 



Just like in the case of the real magnetic field, we assume 
that this is the value of the spin-orbit vector potential at 
the surface of the cylinder. Using Eq. (7), we can rewrite 
the spin-orbit magnetic field and vector potential at the 
surface Eqs. (11) and (14) as 



he 



m e c _ 
4eAf e\r 



A a {R) 



he 

4eA L 



a x n 



m e c 



Tfi 



(15a) 



(15b) 
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respectively. Next, we need to deduce their behavior as 
we move away from the surface towards the interior. 

We consider a cylindrical geometry throughout. We 
expect the spin current to flow within a London pene- 
tration depth of the surface, just like the charge current, 
with the electron spin parallel (antiparallel) to the cylin- 
der axis and antiparallel (parallel) to its orbital angular 
momentum. The 'spontaneous' current of the spin-u car- 
riers is 



Mr) 



(16) 



The total velocity of a carrier of spin a in the presence 
of an external magnetic field is Va,tot = v a + v s . At the 
surface (r = R) the spin current velocity v a (R) is given 
by Eq. (13), and the spin-orbit vector potential by Eq. 
(14) (or (15b)). The equation relating the spin current 
to the spin vector potential is, in analogy with Eq. (5) 



Mr) 



Mr) 



and from Eqs. (16) and (17) 
v a {r) = — 



(17) 



(18) 



analogous to Eq. (8). 

As discussed in ref. Q, the spin current exists due to 
the presence of an electric field E in the interior of the 
superconductor, that satisfies the differential equation[4] 



V 2 (E~E ) 



1 



(E - E ), 



(19) 



where Eo(f) is the electric field arising from the uniform 
positive charge distribution po (V • Eq — Aitpo). E is 
maximum at a distance Ax from the surface and decays 
to zero at the surface (for a cylindrical or spherical sam- 
ple). As discussed by Aharonov and CasherjljJ and other 
authors [l^, the spin-orbit interaction term in the non- 
relativistic limit of Dirac's equation for a particle with 
magnetic moment p, in an electric field can be written 
in terms of a spin-orbit vector potential proportional to 
fl x E. We expect our spin-orbit vector potential to be 
related to the electric field in a similar fashion, and pos- 
tulate that the vector potential driving the spin current 
in Eq. (17) is determined by the difference between the 
electric field near the surface E and the electric field in 
the deep interior Eq through the relation: 

, 1 1 1 rift 

A a {r) = x (E(r) - E (f)) (20) 

eE m 

with E m related to the charge density near the surface 
p_ by the charge conservation condition Eq. (3). 

Eq. (20) reduces to Eq. (15b) for the spin-orbit vector 
potential at the surface, as is seen by setting E(R) = 



and Eq(R) = E m n. From Eq. (18), the spin current 
velocity as function of position is given by 



( f ) = ^x (E(r) 



(21) 



and properly reduces to v® (Eq. (13)) at the surface, as 
required. The spin-orbit vector potential and the spin 
current decay to zero as we move from the surface to- 
wards the interior of the superconductor a distance Al, 
just like the corresponding charge quantities. Taking the 
curl of Eq. (20) (V x (a x E) = a(V ■ E)) and using 
Gauss' law yields 



B a (f) = V x A a (r) 



eE n 



-An(p(f) - p )a (22) 



At the surface r — R 7 p(R) — p-, and if R >> Al wc 
can neglect po in Eq- (22) and obtain: 



BAR) = -^4*,_* 

eE m 



(23) 



which reduces to the un-screened spin-orbit field Eq. 
(15a) upon replacing p_ by its value given by Eq. (3). 
The result Eq. (23) followed from using Gauss' law and 
imposing the global charge conservation condition im- 
plicit in Eq. (3): that the positive charge expelled from 
the interior that gives rise to the maximum electric field 
E m at distance Al from the surface is precisely the same 
as the extra negative charge p_ residing in the surface 
layer of thickness A^. Thus the coincidence of Eq. (23) 
with Eq. (15a) was not built in and supports the consis- 
tency and validity of our framework. 



III. THE FOURTH COMPONENT 

In ref.JH we showed that the electrodynamics in the 
charge sector, under the assumption that the magnetic 
vector potential obeys the Lorentz gauge, could natu- 
rally be formulated in a 4-dimcnsional covariant form in 
terms of the 4-dimensional charge current and vector po- 
tential J = (J.icp), A — (A,icj)) ((j> = electric potential, 
p = charge density). The four-dimensional divergence 
of J and of A vanish due to charge conservation and 
the Lorentz gauge condition respectively. Similarly we 
seek here the fourth component of a 4-dimensional spin 
vector potential A a = {A a ,%<$> a ) by demanding that the 
four-dimensional divergence vanishes: 



DivA 
Using the relation 



1 d 

V-A\ + —ZL<l>a =0. 

c at 



(24) 



V • (a x E) = -a ■ (V x E) 
and Faraday's law we obtain from Eqs. (20) and (24) 



m e cv 



eE n 



(25) 
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and correspondingly we add the fourth component p a 
to form the four-current for carriers of spin a, J a = 
(J a ,icp a ), satisfying the equation (cf. Eq. (17)) 



Pa 



87reS TO A|" 



This equation can be rewritten as 



Cns V °-a ■ B 



2E n 



or as 



4A L £ r 



-ft- B 



(26a) 



(26b) 



(26c) 



and implies that in the presence of an applied magnetic 
field there is excess negative charge corresponding to spin 
direction antiparallell to the applied field B. This cor- 
responds to the electrons near the surface that increase 
their velocity when B is applied. Correspondingly, the 
charge density of the spin component that decreases its 
velocity when B is applied decreases. Note that J a and 
p a are related by the continuity equation DivJ a = 0: 



Ja+ dt 



o 



(27) 



When an external magnetic field is applied, J a near the 
surface acquires a divergence through the induced Vx£, 
which creates the spin imbalance through Eq. (27). 



IV. MAGNITUDE OF EXPELLED NEGATIVE 
CHARGE AND INTERNAL ELECTRIC FIELD 

An applied magnetic field generates a charge current 
near the surface, and according to Eq. (26) it also gener- 
ates a charge density. It is natural to conclude that this 
charge density arises due to the changed carrier velocity. 
From Eqs. (8) and (10), the magnetic field expressed in 
terms of the supcrfluid charge velocity at the surface is 



B = ——v s (R) x n 
and the induced charge density Eq. (26a) is 



,m e c. 

Pa = -(-T-) 

eX L 



8nE, 



-a ■ (v s (R) x h) 



(28) 



(29) 



We argue that it is natural to conclude that the expelled 
charge density p_ will obey the same relation with re- 
spect to the spin current velocity v a that p a bears to the 
charge current velocity v s (Eq. (29)), with an extra fac- 
tor of 2 because of the two spin components contributing 
to the spin current. Hence we postulate that 



or, substituting for p_ using Eq. (3) 

771 _ f ^gf „,0\2 ] 



with solution 



mec„ 
m_ e\ L V ° 



(31) 



(32) 



or, substituting for the spin current velocity from Eq. 
(13) 

he (f>o 

Em = -^>? L = ^W L =2 ^ B (33) 

with 4>q — hc/2\e\ the flux quantum and pb the Bohr 
magneton. Hence the magnitude of the internal electric 
field near the surface is the same as that of the spin-orbit 
effective magnetic field (Eq. (11) or (15a)). 

The value of the expelled charge density near the sur- 
face, p-, follows from Eqs. (3) and (32): 



n s e- 



which can also be written as, using Eq. (13) 



P- = 



n s eh 



4A L m e c 2A L 



or, using Eq. (7), as 



he 
16ne\l 



(34) 



(35) 



(36) 



The charge density induced by the external magnetic 
field can now be written, substituting in Eq. (26a) E m 
by its value Eq. (32) 



Pa 



1 



• B 



8ir\ L 

or, substituting for E m in Eq. 26(b) 

n s e^ 

pa = -T-o- ■ (v s x n) 
zc 



(37a) 



(37b) 



as one would expect by analogy with Eq. (34) . Thus the 
total charge density for spin a is simply 

n s e v°+av s n s eav^ tot 
Pa , ot = _(___) = _ _ (38) 

with iv. tot = v s + cry®, a = +/ — 1 for spin paral- 
lel/antiparallel to the applied magnetic field, and 



Pa=l,tot + Pa=-l,tot 



(39) 



independent of whether or not an external magnetic field 
is applied. 

Equation (5) for the charge current can be written in 
terms of the expelled charge density p_ using Eq. (3) as 



= ( meC )2 V a 



S ■ (i£ x h) 



(30) 



J = P-C 



B x n 



E„ 



(40) 
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where we have used Eq. (10) for the vector potential. 
Eq. (40) has the following interpretation: rather than 
the entire superfluid moving with speed v s (Eq. (4)), we 
may think of the charge current as being carried solely 
by the excess negative charge density p_, moving with 
velocity 



= cB x h/E n 



(41) 



(Note that the total charge density /?_ is not changed 
by the applied magnetic field since p a = —p- a ). The 
speed Eq. (41) would exceed the speed of light if the 
applied magnetic field would exceed E m . However, the 
value we deduced for E m Eq. (33) is just right to prevent 
this from happening, since it is of the order of the lower 
critical field for a type II superconductor [Hj]: 



H c \ — 



47TA! 



(42) 



with k the Ginzburg Landau parameter. Furthermore, 
we argued in Ref.[|[ that superconductivity is destroyed 
when one of the components of the spin current is stopped 
by the applied magnetic field, which corresponds to a 
magnetic field of magnitude E rn . Thus, Eq. (41) can 
be understood as meaning that superconductivity is de- 
stroyed when the speed of the excess charge carriers near 
the surface (p_) reaches the speed of light. 

An alternative and perhaps even more remarkable in- 
terpretation follows from writing the charge current Eq. 
(4) solely in terms of the charge density induced by the 
magnetic field Eq. (37), as: 



J = en s v s = c{p a =-\ - p a =i)v s 



(43) 



Eq. (43) can be read as meaning that the charge 
current induced by the magnetic field is carried solely 
by the induced charge carriers moving at the speed of 
light. Within this interpretation, superconductivity is 
destroyed when the induced charge density of the carriers 
that are slowed down by the applied magnetic field (car- 
riers with spin parallel to the magnetic field) completely 
depletes the component of p~ with that spin orientation. 

In a type I superconductor, superconductivity is de- 
stroyed when the applied magnetic field reaches the ther- 
modynamic critical field H Cl which is smaller than Eq. 
(33) (H c ~ H cl x A L /£ , with £ the Pippard/BCS co- 
herence length). However, in that case the response to 
an applied magnetic field is non-local since the magnetic 
field varies strongly over a coherence length, so these con- 
siderations need to be modified. We will argue in a later 
section that for a type I superconductor E m — H c rather 
than Eq. (33). The expelled charge density p_ is still 
given by Eq. (3), which implies that the interpretation 
of Eq. (41) is the same as discussed above for type II 
superconductors. Fig. 1 shows schematically the spin 
current and charge distribution in a cross section of a 
cylindrical sample. 

The electrostatic energy cost per unit volume for a 
surface charge density p_ expelled from a volume V is 




+ + 4- + 4-' 




FIG. 1: Charge distribution and spin current in a cross section 
of a cylindrical sample. The charge density is negative within 
Xl of the surface and positive in the interior. The spin current 
circulates where the charge density is negative. At the surface, 
the velocity is given by Eq. (13). The horizontal arrows 
indicate the direction of motion of the electrons with magnetic 
moment in the direction indicated by the vertical arrows. 



jiven by (since po is uniform in the interior) 

E 2 



1 

V 



8tt 



-dv 



\3l 

2 8tt 



for the case of a cylinder and 



5 8tt 



(44) 



(45) 



for a sphere, so the electrostatic energy density is smaller 
than the magnetic energy density associated with the 
thermodynamic critical field H 2 /8n for both type I and 
type II superconductors. 

It is interesting to note that Eq. (34) together with 
the charge neutrality condition Eq. (3) and Eq. (7) lead 
to 



m e c 



(46) 



of the same form as the relation between charge velocity 
and magnetic field (Eq. (28)). The kinetic energy density 
of the spin current is, from Eqs. (46) and (7) 



1 p2 

/ 0\2 

-m e (v") n s = — — 

2 K a! 8tt 



(47) 



hence it equals the electrostatic energy density near the 
surface. As is well known, the same relation exists be- 
tween the kinetic energy density associated with the 
charge current and the magnetic energy density near the 
surface. Hence the total kinetic energy density in the 
presence of an applied B and the associated charge cur- 
rent is 



1 



m e (v 2 a 



E 2 . 



B 2 



(48) 
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The charge current screens the external magnetic field so 
it does not get into the superconductor. The spin current 
screens the internal electric field so it does not leak out 
of the superconductor. 



V. SPIN CURRENT SPATIAL DISTRIBUTION 

Having determined the value of the internal electric 
field E m , we return to the formulation of the equations 
governing the spin current. We will slightly generalize 
the equations in Sect. II. As the equation relating spin 
current to spin vector potential we take instead of Eq. 
(17) 



J a {r) ~ JaO = - 



{A„{r) - A a0 (f}) (49) 



with 



and 



A.(r) = X L a x E(r) 



A„ {r) = X L a x E (r) 



•JaO — ~^~ V vO 



(50a) 
(50b) 

(51a) 
(51b) 



The coefficient Al in Eq. (50) results from replacing in 
Eq. (20) the value found for E m , Eq. (32). Note the 
similarity between Eq. (50a) and Eq. (10), which was 
not "built in" . Aside from an inconsequential redefinition 
of A a {A a in Eq. (17) is (A a - A M ) in Eq. (49)), Eqs. 
(49)-(51) differ from Eqs. (17) and (20) in that we allow 
for the possibility of a constant current J CT o deep in the 
interior of the superconductor, where E = Eq. From Eq. 
(49) we obtain using Eqs. (50) and (51): 



and taking the curl 



u x (E(r) - E (r)) (52) 



V x (v a (r) v aQ ) = -J- PW Po J. (53) 



en s X L 
Eq. (53) is approximately satisfied by 



v a (r) = —p(r)a x f 

en s 



(54a) 



The second term in Eq. (55) is zero in the deep interior 
and is smaller than the first term by a factor X^/R near 
the surface, and the first term in Eq. (55) is approxi- 
mately the same as Eq. (53). 

Eq. (54) generalizes the relation found between charge 
density and velocity of the spin current carriers near the 
surface, Eq. (37), to the entire volume. The form Eq. 
(54) is only valid in the absence of applied magnetic 
field, when only the pure spin current exists (w<j = +i = 
—v a= -i). In the presence of both charge and spin cur- 
rent, we write instead of Eq. (54) 



v a (r) = — —p a (r)a x f 

en s 



(56) 



where (modifying the convention used in Sect. IV) we 
denote by p a the total charge density for spin a, pro- 
duced by both the pure spin current and any superim- 
posed charge current induced by an external magnetic 
field (this was called p a ^ot in Sect. IV), and v a denotes 
what was called v a ^tot in Sect. IV Eq. (38). 

In Sect. II we had argued that the spin current should 
die down as we move beyond a London penetration depth 
of the surface towards the interior, by analogy with the 
behavior of the charge current, and formulated the equa- 
tions accordingly. They corresponded to taking J a0 = 
in Eq. (49). However, the finding of a general relation 
between charge density and current in Sect. IV led us 
to include here the J CT o term in the theory, representing 
a counterflowing spin current in the interior induced by 
the charge density po, smaller than the spin current near 
the surface by a factor Xl /R. 

It should be noted however that the addition of the 
constant term J a0 leads to a singularity in the vorticity 
in the deep interior: namely, for r « R 



en s r 



(57) 



diverges as r — > 0. We return to this point in Sect. XII. 

The spin current Eq. (51) expressed in terms of the 
carrier's velocity Eq. (56) takes the form 



J a {r) = -cp a (r)a x f 



(58) 



with 



and denotes the total current carried by the carriers of 
spin <j, including the contribution from the spontaneous 
pure spin current and from the charge current induced by 
an applied magnetic field, if any. Eq. (58) is consistent 
with Eq. (43) deduced for the charge current only. 



VI. SPILL-OVER 



v aQ = v a {r « R) = povxr (54b) 

en s 

since using these expressions we obtain 

^ ,-. -. \ c dp(f) _ c p(r) — po _ . . 
Vx(i(r)-i ) = — '^J-a~ — ' K '- H a (55) 



dr 



en s 



The addition of the constant term J a o to the spin cur- 
rent constitutive relation has another interesting conse- 
quence. Consider the equation for the spin current 



J*{r) - JaO = 



8ttX l 



a x [E(r) - E (r)} 



(59) 
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We argued earlier that at the surface E — 0, Eq — E m f, 
and the spin current velocity is given by the universal 
form Eq. (13). However, these arguments were based on 
the constitutive relation Eq. (17), that did not include 
the Jero term. Eq. (59) no longer satisfies this condition. 
How is this inconsistency resolved? 

Consider a cylinder of radius R. The electric field Eq 
at the surface is 



E (R) = 2irRp h 



(60) 



Assuming the electric field E vanishes at the surface, Eq. 
(59) yields 

en s ^ cp Q _ A c 

—Vo(R) - —o *n = —Rpo(a x n ) (61) 
Bringing the second term to the right side, we have 

= ^-Po(i? + 2A L )(axn) (62) 

which is satisfied if we replace R by R + 2\l in Eq. (2a), 
i.e. if the 'effective' radius of the cylinder is assumed to 
be 



R eff = R + 2\ L 



(63) 



instead of R, which indicates that the negative super- 
fluid 'spills over' a distance 2Al beyond the radius of the 
cylinder. 

To obtain the magnitude of the charge that spilled 
over, note that at r = R the electric field E{R) will no 
longer be zero. To satisfy Eq. (56) with the spin current 
velocity Eq. (13) and Eq = E m (Eq. (33)) requires an 
electric field at R pointing outward, of magnitude 

E(R) = 47rA L p (64) 

which corresponds to a sheet of surface charge density 

VspM = -AlPo (65) 

having spilled out beyond the radius R. Being spread out 
over a radial distance 2Al beyond the surface, it corre- 
sponds to a spill-over volume charge density a sp ui / '(2Xl), 
i.e. 



P spill 



Po 



(66) 



Note also that because of this spill-over effect, the ac- 
tual negative charge density near the surface of the su- 
perconductor is not p-, rather it is p_ + po. Hence upon 
taking the curl of Eq. (20) we need to use 

V • (E - E )r=R = M(p- + Po) - Po) = 47TP- (67) 

This resolves the small discrepancy that we had found 
in Sect. II between the expressions for the unscreened 
spin-orbit magnetic field Eqs. (22) and Eq. (15a). 



Note also that this scenario is consistent with the mi- 
croscopic picture put forth in ref. 8], that electrons ex- 
pand their wavefunctions effective radius from a lattice 
spacing to circular orbits of radius 2Al when a metal goes 
superconducting. For pairs that are right at the surface in 
the normal state, these orbits will expand to a distance 
2Al beyond the body's surface in the superconducting 
state. 

For a sphere of radius R, the total spilled charge is 
4irR 2 a S pui, which using Eq. (2b) and (34) corresponds 
to a number of spilled electrons 



She 



167re 2 Ai 



R 



(68) 



For example, for A^ = 400 A and R — lem this yields 
N S pui = 2 x 10 6 . We had already predicted earlier such 
a spill-over of negative charge beyond the surface of su- 
perconductors based on quite different arguments Q. 



VII. SPIN CURRENT ELECTRODYNAMICS 

We can now simply extend these equations to describe 
the full space and time dependence of the currents and 
charge densities. In terms of the four-vector for spin- 
component a the total current with or without an applied 
magnetic field is given by 



Jcr(r,t) - J aC 



(A a (r,t) - A a0 (r)) (69) 



with 



and 



with 



Mr, t) = (J a (r, t), icp a {r, t)) (70a) 



■ho = (J*o,icp<ra) (70b) 



A a {f,t) = {A a (r,t),ic/) a (f,t)) (71a) 



A M {f) = (A CTO (r),^ CT o(r)) (71b) 



A(j {r, t) = \ L a x E(r, t) + A{r, t) (72a) 



Acro(r) = A L <? x E a (r) 



(72b) 



and 



<Mf,i) = -\ L a-B(r,t) +(/>(r,t) (73a) 



<t><rt{r) = <t>a(r) 



(73b) 



Here, A and <f> are the magnetic vector potential and 
electric potential. Eg and <j>o are the electrostatic field 
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and potential for a uniform charge density po throughout 
the material. The spatial component of Eq. (j6"9")l is 



□ (Jo- — Jo-o) — T2"(<^(T 



>a0 



(80b) 



J<r(r, t) - J a0 = - 



(X L a x E(r,t) + A(r,t)) (74a) 



and the fourth component is 
1 



p<j{r,t) - p a0 



87rA, 



-a-B(f,t) 



8ttA| 



(0(f,t)-0 o (r)) 
( ?4b ) 

The continuity equation sets the four-dimensional diver- 
gence of the four- vector J a equal to zero 



DivJ a = 



(75) 



with the fourth component of the divergence operator 
given by d/d(ict). The Lorentz gauge condition 



C Ot 



(76) 



with 



1 r) 

□ 2 = v 2 - 



2 

c 2 di? 



(81) 



and J a is given in terms of A a by Eq. (69). The equations 
for the charge sector only are simply obtained by defining 
the charge four-current and charge four-potential 



Jc — J<y=+l + J<j=- 



(82a) 



A c = (^ CT=+1 +A,=-i)/2 (82b) 
and similarly for J c q and A c q, and they obey 



J c (f,t) - J c0 = --4r(A c (r,t) - A c0 (f)) (83) 
47rAr 



together with Faraday's law V x E = (—l/c)dB/dt en- 
sures that the four-divergence of A a vanishes, consistent 
with Eq. (69). That A a is a four- vector can be seen from 
the fact that it can be written as 



{A a ) a = —^-€af3jS^pF-yS + A a 



(77) 



with F 7 s the electromagnetic field tensor, A a the 
usual electromagnetic four-vector potential obeying the 
Lorentz gauge condition, and e a p 7 s = +1 (—1) for even 
(odd) permutations of 1234 and zero otherwise. 

The current 4- vectors are given in terms of the velocity 
of the superfluid charge density per spin en s /2, the veloc- 
ity for each spin component v ai and the (excess) charge 
density p a as 

Jcr(r,t) = (^Y-v a (r : t),icp,j(r,t)) (78a) 



J vQ = ( — VaO^CPaQ) 



(78b) 



with v a o given by Eq. (51b) and p a o = Po/2- Using the 
relation (56), they can be written in the remarkable form 



J<j(r,t) = p cr (f,t)c(-a x r,i) 



J*o = x r ' z ) 



(79a) 



(79b) 



which says that the supercurrent density (charge or spin) 
at any point in the superconductor can be understood as 
arising from the excess local charge density moving at 
the speed of light. 

The differential equations determining the behavior of 
all quantities are 



a 2 (A 



A 



1 



<t0) 



T {A a -A M ) 



(80a) 



and Eq. (83), which of course coincide with the equations 
derived in ref.0] for the charge sector. 



VIII. SPIN CURRENT ELECTROSTATICS 

In the absence of time-dependence, the electric field 
in the interior of the superconductor is described by the 
differential equation Eq. (19), and the charge density 
obeys the differential equation Eq. (1). They are related 
by the equation 



E(f)-E (f) =4^A|Vp(r) 



(84) 



For a cylinder, the solution to these equations is given in 
terms of modified Bessel functions I n (z) (Bessel functions 
of imaginary argument) asfloT| 



p{r) = p (l - . ,. p/ . J 
2A L h(iR \ L ) 



E(r) = 



R h{ir/\ L ) 
r h{iR/X L ) 



for a sphere 



1 R 3 sinh(r/X L ) 
p(r) = p (l - - — — TTFTTr-r) 



E(r) 



3A|r f(R/X L ) 



4 R 3 f(r/X L ) 1 

-np r 1 

3 H f{R/X L ) 



with f{x) — xcoshx — sinhx, and for a plane 

R cosh{r / Xl) 



p(r) = p (l 



Xl sinh(R/ Xi 



(85a) 



(85b) 



(86a) 



(86b) 



(87a) 
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E(r) = Airp r[l - 



R sinh(r / Xl) 
r sinh(R/ Xl) 



(87b) 



For R >> Xl they take the following forms: for the 
cylinder 



p(r) = Po(l 



E(r) = 2Tr Po r[l 



R 3 / 2 
2X L r 1 / 2 ' 



(88a) 



R^ 2 

r 3/2 ' 



and Eq = 2irpQr; for the sphere 



p{r) = p (l - 



R 2 
3A L r f 



-(R-r)/X L] 



E{r) = \-K PQ r[l - ^e-^)/^] 
3 r z 



and Eq = (An / 3) pof, and for the plane 



p{r) = p (l - T-e 



-(fl-r)/A I , 



1? 



E(r) = 4np r[l e - {R - r)/XL 



(89a) 



(89b) 



(90a) 



(90b) 



and Eq — Anpof. 

The spin current is given in terms of the electric field 

by 



Mr) 



8irX 



-ax(E- Eq) + Jm 



with 



(91) 



(92) 



JaO = ~2 P °° X ^ 

For the cylinder, Eq. (91) yields 



= E Y [ ^ e ~ lR ~ r)/XL - l]d x f (93) 

and a similar expression for the sphere. Both can be 
written as 

Mr) = -|[p-e~ (fl - r)/At + po]ff x f = (94) 

or 

J a (r) - J a (R)e-( R - r V XL + J a0 (95) 
The magnitude of the spin current near the surface is 



*Amp , 400A 

1° v ( 

cm 2 A L (1) 



•i 3 



J a (R) = 1.02 x 10 s f x ( — ) (96) 



and in the deep interior 



2X L 
~R 



MR). 



(97) 



The direction of the spin vector a is arbitrary in the 
spherical geometry. If an external magnetic field is ap- 
plied in the z direction, a charge current develops ac- 
cording to Eq. (5) in the azimuthal direction and defines 
the spin quantization axis. Then the up and down spin 
components of the total current are simply obtained by 
using the z axis as the spin quantization axis and adding 
the charge and spin current contributions. Similarly, a 
magnetic field applied parallel to a planar surface will 
polarize the spin current with the spin quantization axis 
parallel to the applied field. 

If an external electric field E ex t is applied, Eq. (19) 
remains valid but the boundary condition changes [14j. In 
a spherical or cylindrical geometry there is no external 
electric field in the absence of applied electric field, hence 
we have at the surface 



E = E e . 



(98) 



Assuming there are no 'normal electrons', the electric 
field is continuous at the surface and penetrates a dis- 
tance Al[4|. So right at the surface we have 



Ja(r) = -3— x (E ext - Eq) + J ( 



8irX 



(TO 



(99) 



Consequently, an external electric field that points in 
(antiparallel to Eq) will enhance the spin current, and 
one that points out will suppress it. If the magnitude of 
the electric field applied is larger than the critical field 
E m (eq. 33) it will completely suppress the spin current 
and hence destroy superconductivity. 

If an electric current circulates along the axial direction 
of the cylinder, one expects that in the absence of applied 
electric field that above a critical current 



_ c H c l 

J cl — -. r — 

4tt X l 



(100) 



vortices start to penetrate the superconductor [13]. If an 
electric field normal to the surface is applied, this expres- 
sion will be modified to 



Jcl — 



c H cl ± E ext 



4tt 



Xj 



(101) 



with the + (— ) sign corresponding to electric field point- 
ing into (out of) the cylinder. Thus, for a cylinder placed 
between the plates of a capacitor the flux penetration for 
a given current flowing along the cylinder should be af- 
fected by the electric field of the capacitor and depend on 
the position in the surface relative to the capacitor plates: 
the flux penetration should be larger for the region near 
the negative capacitor plate. 
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• <B„> 





points radially outward, 



(a) Discrete case 



(b) Continuous case 



FIG. 2: (a) The equivalent electric dipole for an electron with 
magnetic moment pointing out of the paper propagating up- 
ward is attracted by the ion on the left and repelled by the 
ion on the right, resulting in a force pointing left, (b) The 
torque from the radially outward pointing electric field of the 
ions changes the angular momentum of the electron so that 
it deflects to the left, aquiring angular momentum pointing 
out of the paper, i.e. parallel to p. Similarly an electron 
with magnetic moment pointing into the paper deflects in the 
opposite direction. 



IX. ORIGIN OF THE SPIN CURRENT 

As discussed in Ref. Q , the physics that we are propos- 
ing as underlying both the generation of spin current and 
the Meissner effect in superconductors is an expansion 
of the electronic wavefunctionQ, from a lattice spacing 
to a radius 2Al, as the system enters the superconduct- 
ing state. In the presence of a magnetic field, the az- 
imuthal Lorentz force on radially outgoing electrons gives 
rise to the charge current that 'expels' the magnetic field. 
An azimuthal spin-orbit force originating in the positive 
background 15] gives rise to the spin current, whether or 
not an external magnetic field is present. 

The fact that a magnetic moment propagating in a 
periodic lattice experiences a net transverse force was 
proposed in Ref. [151 ] as an explanation for the anoma- 
lous Hall effect in ferromagnetic metals. An electron with 
magnetic moment p, propagating with velocity v is equiv- 
alent to an electric dipole [If 



P 



7- X p 

c 



(102) 



with 7 = (1 — v 2 /c 2 ) 1 ~ 1 for non-relativistic speeds. 
An electric dipole in an electric field experiences a force 
F = V(p ■ E) and a torque f = p x E. As shown in Fig. 
2(a), for a discrete positive charge distribution the force 
acting when the electron passes through the line joining 
two ions is in direction jxxv. Averaging over the unit cell 
we showed in Ref. [ll| that the net force is proportional 
to the averaged second derivative of the lattice potential 
and in direction p x v, corresponding to an 'effective' 
magnetic field in direction parallel to p. 

For the 'mesoscopic' scale under consideration (2\ L ) 
it is more appropriate to consider a uniform, positive 
charge density, as shown in Fig. 2(b). The electric field 



E = 2irpr 



(103) 



with p the ionic charge density. The torque by the electric 
field acting on the equivalent electric dipole changes the 
angular momentum 

dL d ._, _« _. ,v . => . . 

— = m e — (r x v) = t = (- x p) x E (104a) 
dt dt c 



hence 



dv 
~dt 



-2-irprx (- x p) 



(104b) 



which is equivalent to the effect of an 'effective' magnetic 
field in direction parallel to p 



-2n—p 

e 



(105) 



on the electron, also causing it to deflect in direction 
p x v, as shown in Fig. 2(b). Eq. (105) is the same as 
Eq. (11), since the superfluid of negative charge density 
en s should respond to the positive ionic charge density 
p = -en s . 

The corresponding 'unscreened' vector potential, de- 
fined by V x A a = Bfj is given by 



- Bfj X T p _ ^ 

A a = — - — = -ir-p x r 
2 e 



which yields 



Am„c 



-a x E. 



(106) 



(107) 



for E given by Eq. (103). The spin-orbit interaction term 
in the non-relativistic limit of the Dirac equation for an 
electron in an electric field E is 

H s .o. = -— %-s3> (E xp) = — p-A a (108a) 
4mgC m e c 

where p is the momentum operator and A a is given by 
Eq. (107). Thus, we can understand the relation Eq. 
(18) between velocity and spin-orbit vector potential as 
arising from minimization of the energy-momentum re- 
lation that arises from the non-relativistic limit of the 
Dirac equation 



E(p) 



P 

2m e 



-P- A a 



(108b) 



with p — m e v a and A a given by Eq. (107). 

It should be mentioned however that our conclusion is 
opposite to that of Chudnovsky[l9| concerning the sign 
of the spin-orbit force acting on an electron moving in a 
positive background. Ref. [191 ] defines a spin-orbit vector 
potential as Eq. (107) with opposite sign, and concludes 
that the force is in opposite direction to that depicted 
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in Fig. 2. We argue that it is physically clear that the 
electron will deflect as shown in Fig. 2, so that the elec- 
tron orbital angular momentum will point antiparallel to 
its spin angular momentum to minimize the spin-orbit 
energy, just like in atoms. The difference between our 
conclusion and that of ref. [l!| can be accounted for by 
the momentum of the electromagnetic field, which was 
calculated by Aharonov et al as [20| 

Pfieid = ^ I E x B = -E x /x (109) 

and points in opposite direction to the mechanical mo- 
mentum of the electron in Fig. (2) and is twice as large. 

As electrons of opposite spin move radially outward a 
distance 2A^ they acquire azimuthal velocity 



A a (r = 2A L ) 

m e c 



-a x r 



4m e A L 

where we have used p = |e|n s [8( and Eq. (7). 



(13) 



electrons per unit volume in such orbits, the magnetic 
susceptibility per unit volume is 



X 



d 2 AE 



dH 2 Am e c 2 
and for orbits of radius r = 2Al, as proposed here, 



fill) 



X 



n s e 2 
2 l 



1 

47T 



(112) 



using Eq. (7) for the London penetration depth. Eq. 
(112) is the condition for perfect diamagnetism. In- 
stead, the Landau diamagnetic susceptibility in the nor- 
mal state is given by Eq. (Ill) with r = kp, with fcp 
the Fermi wavevector. 

Similarly, Eq. (110) yields that the increase in energy 
per unit volume for n s electrons per unit volume in orbits 
of radius r — 2Al is 



u^n s AE=-^(2\ L ) 2 H 2 
Hm e c 



8tt 



(113) 



X. PHYSICAL INTERPRETATION 

There are three key elements underlying the physical 
picture considered here that relate to the theory of hole 
superconductivity (2lj| : (i) The wavefunction in the nor- 
mal state has to be confined to short distances, which 
requires an almost full band (hole transport in the nor- 
mal state); (ii) The carriers forming the superfluid are 
long wavelength e/ectronsQ, with negative charge^^; 
(iii) In quantum mechanics, expansion of the wave func- 
tion is associated with lowering of kinetic energy, which 
is the driving force for hole superconductivity [231] . 

It is interesting to note that an early explanation of su- 
perconductivity proposed by J.C. Slater in 1937(24) had 
some key elements in common with what is being pro- 
posed here. Slater proposed a model where "the wave 
functions correspond to electrons which can wander for 
some distance through the metal" , and showed that "to 
produce superconductivity the orbits must be of order of 
magnitude of 137 atomic diameters". 

Consider an external magnetic field H applied to an 
electron orbiting with radius r in a plane perpendicular 
to H . Classically, in changing the field from H to H+dH 
the electron changes its speed by dv = (e/m e c)(rdH/2) 
due to Faraday's law, and its orbital magnetic moment by 
rf/i = (er/2c)dv antiparallell to the field, resulting in an 
energy increase dE = Hdp,. When the field is increased 
from to H, the total energy increase is 



AE 



H e 2 
o 4m e c 



2 tt! 

r H 



8m e c 



;r 2 H 2 



(110) 



The same result is obtained quantum-mechanically, with 
r 2 replaced by the mean square radius in the plane per- 
pendicular to H, r 2 — >< x 2 + y 2 >. If there are n s 



again using Eq. (7), so that the system will remain in 
this state as long as the 'condensation energy' density of 
the state is greater than the energy cost H 2 /8ir. This is 
of course the condition that defines the thermodynamic 
critical field -EL [131]. 

Slater's arguments[24] proceeded similarly. He pointed 
out that the radius of the orbit in the susceptibility ex- 
pression Eq. (Ill) should be such that \ = — 1/4tt for 
perfect diamagnetism; taking n s = 1/rf 3 , with the vol- 
ume per superconducting electron, and setting d = 2ao, 
with ao = h 2 /m e e 2 the Bohr radius, yields 



d 2 



m e c 2 d 2 ,hc,~ 

-|-z = zbr) 



(114) 



hence Slater concluded that "the orbits must be of order 
of magnitude of 137 (= hc/e 2 ) atomic diameters". 

Furthermore, Slater proposed as a criterion for the crit- 
ical magnetic field that will destroy superconductivity to 
compare the Landau level spacing of energy levels for 
fixed z-component of momentum in an external field H 



AE U = 

m„c 



(115) 



to the spacing of discrete energy levels in a 'box' of size 
137 atomic units. Instead, we argued in ref.Q that su- 
perconductivity will be destroyed when the external mag- 
netic field completely stops the spin current orbital mo- 
tion of the electron with spin antiparallell to H, i.e. for 



H = B n 



(116) 



with B a given by Eq. (15a). The Landau level energy 
spacing for such a magnetic field is 



AE U = hv c 



eh 
m e c 



-B„ 



Am e X 2 L 



(117) 
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which is of the same order of magnitude as the energy 
level spacing of electrons confined to a region of size 
~ 2Al- Thus our criterion is essentially equivalent to 
Slater's criterion. An equivalent form of this criterion is 
that superconductivity will be destroyed when the energy 
increase for electrons in a magnetic field H for orbits of 
radius r = 2Al, Eq. (110), becomes comparable to the 
spacing of energy levels in a box of such size. 

Of course Slater's scenario left several questions unan- 
swered that our present scenario addresses, namely: (i) 
what happens to those orbits when there is no external 
magnetic field applied? As we showed in [8J, orbits with 
radius 2Al arise from the spin-orbit interaction under 
the constraint that the orbit should enclose precisely one 
flux quantum of spin orbit flux, so that the wave func- 
tion is single-valued. And (ii) how do the orbits arise 
when a metal is cooled into the superconducting state, 
and how is the Meissner current generated? It is the pro- 
cess of expansion of the orbits from negligible radius (of 
order of a lattice spacing) to radius 2Al that generates 
dynamically the azimuthal velocities required. 

It is also interesting to note that the physics of super- 
conductivity discussed here provides a natural 'classical' 
explanation for the existence of macroscopic phase coher- 
ence in superconductors. The fact that electrons traverse 
overlapping orbits of radius 2Al implies that the orbits 
of many electrons cross the orbit of any given electron, so 
the motions needs to be synchronized: to avoid collisions 
the angular position of the i-th electron 9i(t) = uj c t + 9io 
has to bear a definite relation with the angular positions 
of the electrons in the overlapping orbits, that persists 
over time. It is also not possible to change one orbit 
without affecting the synchronization with all the others, 
which gives an intuitive understandingto the notion of 
'rigidity' of the superconducting state[5j. Finally, it is 
interesting to note that several workers in the pre-BCS 
era proposed mechanisms of superconductivity based on 
the notion of 'spontaneous currents' that would exist in 
the superconductor in the absence of applied external 
fields [25[. However these were charge currents rather 
than spin currents. 

XI. ANGULAR MOMENTUM 

Within our model each electron in the supcrfluid moves 
in a circular orbit of radius 2Al, with electron spin per- 
pendicular to the plane of the orbit in direction antiparal- 
lel to the orbital angular momentum, and velocity given 
by Eq. (13). Hence the mechanical orbital angular mo- 
mentum of the i-th electron is 

li.orb = m e (2\ L )r x v a = --a (118) 

where the last equality follows from Eq. (13). The elec- 
tron spin angular momentum is 

h ^ 



and consequently its total angular momentum (orbital 
plus spin) is exactly zero. Thus, the transition to su- 
perconductivity can be thought of as a "quenching" of 
the electron's spin angular momentum by the develop- 
ment of an opposite-pointing orbital angular momentum 
of the same magnitude. 

Let us consider the total mechanical angular momen- 
tum for electrons of spin a in a superconducting cylin- 
der of radius R and height h, choosing as usual the spin 
quantization axis parallel to the cylinder axis. The total 
number of electrons with spin a for superfluid density n s 
in the entire volume V = irR 2 h is 

n 

AW = -j v (120) 

Hence the total mechanical angular momentum carried 
by electrons with spin a is the product of Eq. (120) and 
(118): 

— > — * Tha ft 

L V oi,a = N vo i t(T li t0rb = -n s Vm e \ L v a a = 

On the other hand, the orbital angular momentum for an 
electron of spin a in the surface layer of thickness Al is 

-* R -* 

hurf = m e Rf X v a = ——h t0 rb (122) 

and the total number of electrons of spin a in the surface 
layer of thickness Al is 

N surf , a = 2nR\ L h^ = = ^N vol>a (123) 

so that the total mechanical angular momentum carried 
by the electrons of spin a in the surface layer is 

Lsurf,a N sur f a l sur f N vo i a li or b L vo ifj : (124) 

the same as Eq. (121). So we may think of the total 
mechanical angular momentum in two equivalent ways: 
(i) each superfluid electron in the bulk carries mechanical 
angular momentum h/2, or (ii) the superfluid electrons in 
the surface layer of thickness Al each carry a mechanical 
angular momentum RK/AXl (resulting from their orbital 
motion with speed v a and radius R), and the electrons 
in the interior carry no mechanical angular momentum. 
The total mechanical angular momentum of the super- 
fluid electrons of each spin component in the surface layer 
equals in magnitude the total spin angular momentum 
of the superfluid electrons of that spin component in the 
entire volume. 

If we cool the system into the superconducting state in 
the presence of an applied magnetic field B, the electrons 
of both spin orientations acquire an extra contribution to 
their orbital angular momentum as they expand to their 
orbits of radius 2Al- The azimuthal velocity acquired is 

vj, = —f x B (125) 
m e c 
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due to the action of the Lorentz force as the orbits ex- 
pand to radius 2Al[|[. increases the speed of the elec- 
trons with jl parallel to B, and decreases the speed of 
the electrons with pi antiparallel to B. The extra angular 
momentum per electron due to the magnetic field is 

2eA 2 - 

l B = m e (2X L )f x v<p = -^B (126) 

and the total extra angular momentum acquired by the 
electrons is obtained by multiplying I b by the total num- 
ber of electrons: 

c 2e 

as expected, since 

? 2m e c _ , „„. 

L e = — m (128) 

e 

for electrons, with m = nR 2 hM the induced magnetic 
moment and M = —B/4ir the required magnetization 
to cancel the magnetic field in the interior. Hence, just 
like for the angular momentum in the spin current, the 
angular momentum in the Meissner current in the sur- 
face layer can be interpreted as arising from the orbits of 
radius 2A^ of each electron in the bulk. 

How is this electronic angular momentum compen- 
sated? We raised this puzzle in ref. [26]: in the conven- 
tional theory of superconductivity the charge distribution 
is homogeneous, so there is no angular momentum in the 
electromagnetic field. The electrons in the Meissner cur- 
rent carry mechanical angular momentum L e , but there 
is no mechanism in the conventional theory to generate 
a compensating angular momentum of the ions. 

This conundrum is naturally resolved in the present 
scenario. Recall that in the absence of magnetic field 
the 2Al circular orbits arise from the interaction of the 
magnetic moment moving outward with the positive ionic 
lattice. The physics is illustrated schematically in Fig. 3. 
The ions exert a torque on the electron 

T ie =pxEi (129) 

where Ei is the ionic electric field and p is the equivalent 
electric dipole Eq. (102). As the ions exert the torque 
on the electrons the ionic lattice is subject, by Newton's 
third law, to an equal and opposite torque: 

Tei = ~Tie (130) 

In the absence of applied magnetic field, r e j is equal and 
opposite from electrons with spin up and spin down, 
hence there is no net torque on the ions. However, in 
the presence of B the velocity of both spin electrons is 
modified according to Eq. (125), and hence the net f ei is 
no longer zero: as the electrons acquire the extra angular 
momentum Eq. (126), the extra torque on the ions points 
opposite to the direction of B and generates an equal and 




FIG. 3: The thick line shows the trajectory of an electron 
of magnetic moment jl pointing out of the paper starting at 
the origin till it reaches the circumference of radius 2Al. The 
equivalent electric dipole moment p points perpendicular to 
the trajectory and outward, the ionic electric field Ei points 
radially outward and the cross product p x Ei is parallel to 
jl. The torque exerted by the lattice on this electron fi e = 
px Ei points out of the paper and goes to zero as the electron 
reaches the circular trajectory. For the electron of opposite 
spin, the torque points into the paper and the trajectory is 
a mirror image across the vertical direction. In the presence 
of a magnetic field pointing out of the paper the torque on 
this electron becomes larger and the one on the electron of 
opposite spin becomes smaller. Hence there is a net torque 
exerted by the ions on the electrons pointing out of the paper, 
and a reaction torque exerted by the electrons on the ions 
pointing into the paper that causes the body as a whole to 
rotate clockwise. 

opposite angular momentum for the lattice, Uons = —Ib 
per electron, and a total ionic angular momentum 

U = -L e (131) 

so it would appear that the total angular momentum is 
conserved. 

However that is not quite the whole story. Because 
the resulting charge distribution is inhomogeneous, there 
will also be some angular momentum stored in the elec- 
tromagnetic field: 

L fleld = ^ c J d 3 rr x (E x B) (132) 

which we can estimate as 

V\ L 

L field = — -z — E m B (133) 
2-kc 

since the region where both electric and magnetic fields 
are nonzero is only the surface layer of thickness ~ A^. 
Using Eqs. (130), (124), (33) and (7) 

L field = 5 \ L e = -L(^l)L e ~10- 6 £ e (134) 

Sm e CAL 07T \l 
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with A c = h/m e c the Compton wavelength. Hence 
we conclude that the spin-orbit interaction can account 
for 99.9999% of the angular momentum conservation 
puzzlefli] through Eq. (131), but we are still missing 
0.0001%! 

This 'missing' angular momentum is of course the tiny 
bit of extra electronic angular momentum that is ac- 
quired by the electrons in p_ that moved outward cutting 
magnetic field lines near the surface in the process. The 
change in azimuthal velocity of an electron near the sur- 
face moving outward a distance Xl in a magnetic field B 
due to the Lorcntz force is 



Av = AA 

m e c 



A(f> eX L B 



2nR 



m e c 



and the acquired angular momentum is 



I = — RX L B 

c 



(135) 



(136) 



parallel to the magnetic field. The number of electrons 
that moved out to the region within A^ of the surface 
and in the process cut through magnetic field lines can 
be estimated as 



N e ,out = —2nRX L h 
e 



(137) 



hence the mechanical angular momentum gained by the 
outflowing p_ charge is 



L e — N eout l — 



2p-VX\ § _ 



VXj 
~2ttc 



-E m B 



(138) 



which is equal and opposite to the angular momentum 
of the field Eq. (133). (Note that Eq. (2a) was used in 
deriving Eq. (138)). Consequently, angular momentum 
is conserved when a metal is cooled into the supercon- 
ducting state in the presence of a magnetic field, since 
the total angular momentum above T c is zero, and from 
Eqs. (131), (133) and (138) 



L e + L e + L fi e id + Li — 0. 



(139) 




(a) Type I superconductor 




(b) Type II superconductor 




(c) Cross-over 



FIG. 4: Type I versus type II materials. £ is the distance be- 
tween the centers of the orbits of magnetic moment up elec- 
tron (left circle) and magnetic moment down electron (right 
circle). The small dots denote the centers of the orbits, the up 
and down arrows the direction of the electron's magnetic mo- 
ment, their direction of motion is indicated with a horizontal 
arrow. In type I materials ((a)), £ > 2\l and the orbits don't 
overlap. In type II materials ((b)) with £ > 2A_t, a normal 
vortex core can be enclosed by both orbits, (c) denotes the 
cross-over situation, with £ = 2Xl 



XII. TYPE I VERSUS TYPE II MATERIALS 

In a type II superconductor, magnetic flux penetrates 
the body for H c > H c \ and divides itself into filaments, 
each carrying a flux quantum fa. In the vortex core 
of size £, the 'coherence length', the system is normal. 
Instead, in a type I superconductor the flux is excluded 
until the applied field reaches the value H = H c , the 
thermodynamic critical field, and for H > H c the entire 
material becomes normal. The thermodynamic critical 
field is given by 

/3 he Xl 
9^nr~ Hcl T (140) 

with £o the Pippard-BCS coherence length [l3l|. A sim- 
ilar relation holds at finite temperatures between the 



Ginzburg Landau coherence length £ and the penetration 
depth. We will not distinguish here between £ and £o- For 
type II superconductors, £ < Xl, hence H cl < H c . 

We propose the following interpretation of type I ver- 
sus type II behavior, depicted in Fig. 4. Up and down 
spin electrons have orbits of radius 2A_l, but the orbits co- 
incide only in the extreme type II limit where Xl >> £■ 
£ represents the (average) distance between the centers 
of the orbits of spin up and spin down components of the 
Cooper pair. 

In a type I superconductor (Fig. 4(a)), £ >> Xl and 
the orbits of ] and j. electrons are disjoint. Because or- 
bits are time-reversed partners, a magnetic field cannot 
thread one of the orbits and not the other. Instead, in a 
type II material with £ << Xl (Fig. 4(b)) magnetic flux 
can be enclosed by both orbits simultaneously. When 
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(a) Type I 




and obtain 



(b) Type II 



FIG. 5: Type I versus type II materials in a magnetic field 
(schematic), (a) The wavy lines connect the centers of the two 
members of a Cooper pair. In the intermediate state, normal 
regions exist with laminar or other shapes (shaded regions) of 
size larger than £ where Cooper pairs are destroyed and the 
magnetic field is H c . (b) £ << Xl and the orbits of the two 
members of a Cooper pair almost overlap. The vortex cores 
(shaded regions) are within the orbits of the two members of a 
Cooper pair and have magnetic field H c \ when vortices don't 
overlap. 



H 



he 



(142) 



which is essentially the thermodynamic critical field Eq. 
(140). Thus, in a type I material superconductivity is 
not destroyed when the magnetic field stops the spin 
current Q (since H c i > H c here) but rather when the 
diamagnetic energy cost is large enough that the pertur- 
bation significantly mixes the unperturbed energy levels, 
i.e. disturbs the wave function so that it is no longer 
'rigid'. 

What is the expelled charge density p_ in type I ma- 
terials? The conclusion that E m is given by Eq. (33) 
(E m ~ H c \) is not valid because the electrodynamic re- 
sponse is non-local. Furthermore, the electrostatic en- 
ergy cost Eq. (44) would be much larger than the con- 
densation energy in this case. The amount of p_ 'needed' 
to sustain the Meissner current up to H = H c is, accord- 
ing to Eq. (41), simply p_ = —H c /4:iT, so we conclude 
that for type I materials 



E m — H r 



(143a) 



the magnetic field first enters a type II superconductor, 
each orbit will enclose a single flux quantum 4> and the 
magnetic field at the center of the vortex is H c i . In the 
mixed state, at the core of radius £ = distance between 
the centers of the orbits, the magnetic field can have any 
value between H c \ and H C 2 and there will be several vor- 
tex cores within the orbits of a Cooper pair with the total 
flux threaded by each orbit an integer multiple of 0o- The 
cross-over between type I and type II regimes occurs for 
£ = 2Al (Fig. (4c)), where the overlapping part of the 
orbits can just enclose the vortex core. When the centers 
of the orbits are at distance £ > 2Al (type I), the mag- 
netic flux will destroy the Cooper pairs altogether, but 
depending on the sample shape there can be intertwined 
regions of normal and superconducting phases (interme- 
diate state). 

Fig. 5 shows schematically type I and type II super- 
conductors in a magnetic field. In type I superconductors 
with non-zero demagnetizing factor there will be normal 
regions usually of laminar shape intertwined with the su- 
perconducting regions. The uncompensated orbits in the 
neighborhood of the normal regions and vortices will give 
rise to spin current and excess negative charge similar to 
the behavior near the surface. 

This interpretation also provides a rationale for the 
form of the critical field Eq. (140). For £ > \ L , the 
'confinement' region for an electron in the Cooper pair is 
£ rather than 2Xl- Following Slater's reasoning [24j. we 
equate the diamagnetic energy cost for the electron in an 
orbit of radius 2Al in the presence of an external field H 
to the spacing of energy levels in a box of length £ : 

AE «^^ 2H2 ^ (141) 



Hc ft l he 

p - = -4^- L = \2^^ir (143b) 

rather than Eqs. (33) and (34) applicable to type II 
materials. Because of the non-locality we expect the spin 
current density near the surface in type I materials will 
no longer be given by Eq. (16), but rather 

4 = -yffxr (144) 

In type I materials Eq. (144) is smaller than Eq. (16) by 
a factor A^/£. 

Next we address the fact that disorder is known to in- 
crease the value of Xl and turn a type I into a type II 
material. This can be easily understood in the picture 
discussed here. In the presence of defects that weaken 
or destroy superconductivity locally, the superconductor 
will expel negative charge towards those regions, thus de- 
pleting the magnitude of p_ near the outer surface and 
decreasing the ability of the superconductor to shield ex- 
ternal magnetic fields. These internal weak regions (grain 
boundaries, defects, etc) will have excess negative charge 
and accompanying spin current orbiting around them, as 
shown schematically in Fig. 6. If we cool this supercon- 
ductor in the presence of a pre-existent magnetic field, 
it will be unable to expel the magnetic field from those 
regions, thus resulting in an incomplete Meissner effect. 

This picture also suggests that in extreme type II ma- 
terials there will be an inhomogeneous charge distribu- 
tion and spin current distribution arising from inhomo- 
geneities and disorder also in the absence of applied mag- 
netic field. Such inhomogeneities are observed in the 
underdoped regime of high T c materials f27j. which are 
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FIG. 6: Schematic depiction of a superconductor with strong 
disorder in the absence of applied magnetic field. Defects, 
grain boundaries, vacancies, etc. will result in patches of 
normal regions (hatched areas) surrounded by spin currents 
(dashed lines, with arrows pointing in the direction of flow of 
electrons with magnetic moment pointing out of the paper) 
and excess negative charge density (gray areas). The figure 
also shows the excess negative charge and spin current near 
the surface and the spill-over beyond the surface (denoted 
by the full line). The equivalent electric dipoles (Eq. (102)) 
point outward near the outer surface and towards the normal 
regions in the interior. If the system is cooled in the pres- 
ence of a magnetic field, magnetic flux will be trapped in the 
hatched regions and a charge current will flow around those 
regions together with the depicted spin currents. The small- 
est normal regions have diameter of a coherence length and 
enclose one flux quantum of spin-orbit flux in the absence of 
applied magnetic field. . 



XIII. SUPERCONDUCTING MATERIAL 
PARAMETERS 



From experimentally measured values of the penetra- 
tion depth and critical fields we can infer the magnitude 
of the maximum internal electric field E m , surface charge 
density p_, velocity of the spin current carriers at the 
surface, and the spin current density flowing near the 
surface. We use the expressions: 



H c i — 



he _ 1.64 x 10 8 
Ai(A) 2 



4eA| 



Gauss 



E m = H cl {G) x 300 V/cm 
for type II materials (H c i < H c ) and 

E m = H C {G) x 300 V/cm 
for type I materials. 

m e c 2 



47re 2 A| A L (A) 2 A 3 



(145) 



(146a) 



(146b) 



(147a) 



n s (electrons/ion) = n s 



A 



' 0. 602 p(gr/ cm 3 ) 
with A the average mass number of the compound, 



(147b) 



AttXlc 



= 5.524 x 10 



_ 11 E m (V/cm) electrons 



*l(A) 



A 3 



(148) 



expected to be in the extreme type II limit (largest Xl, 
smallest £) 28] . The existence of spin currents associated 
with regions of charge inhomogeneity predicted by our 
model has not yet been experimentally tested. 

Finally, recall that we argued in Sect. V that the equa- 
tions for the spin current spatial dependence should in- 
clude the constant term J aa , giving rise to a singularity 
in the deep interior. Presumably this implies that even in 
the absence of disorder at least one 'vortex', i.e. a normal 
region of size £ surrounded by spin current, has to exist 
in the deep interior of any type II superconductor. This 
would imply that any type II superconductor is topolog- 
ically a torus of genus larger or equal to 1 even in the 
absence of applied magnetic field, hence that a Meissner 
effect with 100% flux expulsion can never be attained in a 
type II superconductor no matter how small the applied 
external field: at the minimum one trapped magnetic flux 
quantum will always remain in the interior. 



h 2.896 x 10 7 , „„. 

v a = — = — — T — cm s (149) 

4m e A L X L (A) 



J spin J <j— + 1 J(T——1 — P—C 

= — (-^r) x 4.802 x 10 15 Amps/cm 2 (150) 
e A 3 

We list the values for a variety of materials including 
type I and type II superconductors in Table I. From this 
table we learn that materials with the largest values of 
P-, E rn and J spin are Nb, Pb and Hg. They are close 
to the cross-over between type I and type II regimes and 
have relatively high T c . Hence these materials should be 
favored in experiments aimed at detecting this unconven- 
tional physics. Instead, the higher T c materials known to 
date (last entries in Table I and similar) are strongly type 
II and the magnitude of the spin current and associated 
quantities is substantially smaller. 
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TABLE I: Properties of superconducting elements and compounds. The superfluid carrier density n s is extracted from the 
measured penetration depth \l using expression Eq. (7), so it corresponds to what is usually interpreted as n s m e /m* , with 
m* the effective mass. The spin current density J 3p i„ is defined by Eq. (144) and includes the contribution from both spin 
components. 



Material 


T C (K) 


A L (A) 


H C (G) 


H cl {G) 


n a {e/A 3 ) 


n s (e/i) 


10 6 p-/en s 


E m {V/cm) 


v a {cm/s) 


J spin (10 6 A/cm 2 ) 


Cd 


0.56 


1300 


30 


97 


0.00167 


0.0361 


0.23 


9,000 


22,275 


1.83 


Zn 


0.875 


290 


53 


1955 


0.0336 


0.510 


0.090 


15,900 


99,855 


14.6 


Al 


1.14 


500 


105 


658 


0.0113 


0.188 


0.31 


31,500 


57,916 


16.7 


In 


3.40 


640 


293 


401 


0.0069 


0.180 


1.10 


87,900 


45,247 


36.4 


Sn 


3.72 


510 


309 


632 


0.0109 


0.293 


0.92 


92,700 


56,780 


48.0 


Hg 


4.15 


410 


412 


978 


0.0168 


0.412 


0.99 


123,600 


70,629 


80.2 


Pb 


7.19 


390 


803 


1080 


0.0186 


0.561 


1.83 


240,900 


74,251 


163.7 


Nb 


9.50 


400 


1980 


1028 


0.0176 


0.324 


2.41 


308,400 


72,395 


205 


MgB 2 


39.2 


1800 


2600 


51 


0.00087 


0.0084 


0.54 


15,300 


16,088 


2.26 


Lai.BsSr .15CM2O4 


37.3 


2500 


2370 


26 


0.00045 


0.0071 


0.38 


7,800 


11,583 


0.826 


YBa 2 Cu 3 7 


90 


1500 


6450 


73 


0.0013 


0.017 


0.62 


21,900 


19,305 


3.88 



XIV. SUMMARY AND DISCUSSION 

In this paper we have extended the electrodynamic 
equations for the charge sector of superconductors pro- 
posed in Ref.[3] to describe the electrodynamics of the 
spin sector, based on the proposal of ref. [8| that a spon- 
taneous macroscopic spin current flows within a London 
penetration depth of the surface of superconductors, in 
the absence of applied external fields, with carrier speed 
given by the universal form Eq. (13). We have shown 
here that the space and time dependence of the spin cur- 
rent can be described within the same four-dimensional 
framework developed earlier to describe the charge elec- 
trodynamics within the theory of hole superconductivity, 
which predicts that negative charge is expelled from the 
interior to the surface when a metal goes superconduct- 
ing. The extension of the theory to describe the spin 
sector is necessitated by the fact that an interior elec- 
tric field exists in superconductors within this theory. 
Instead, in the conventional London theory there is no 
electric field in the interior of superconductors, hence no 
separate electrodynamic description of the spin sector is 
required. 

The formalism led us to the determination of the mag- 
nitude of the expelled negative charge, which had been 
left undetermined in our earlier work. Remarkably, the 
maximum electric field in the superconductor resulting 
from the charge expulsion was found to have identical 
magnitude as the spin-orbit field that gives rise to the 
spontaneous spin current, and to have similar magnitude 
as H c i, the conventional lower critical field of type II 
superconductors. For type I superconductors, we con- 
jectured that the maximum electric field in the interior 
is H c . In addition, we found a remarkable connection 
between the magnitude of excess charge at any point in 
space and time and the speed of the superfluid at that 
point, 

en s 

PaC=—V a (151) 
which says that any superfluid current (charge or spin) 



can be understood as originating in the excess charge 
moving at the speed of light. We also found a remarkable 
parallel between the kinetic energy of the carriers due 
to charge and spin currents and the respective energy 
densities of magnetic and electric fields near the surface: 

\m e {v s ) 2 n s = — (152a) 



-m e (v°fn s = (152b) 

All of this hints at deeper underlying physics that re- 
mains to be uncovered. 

The macroscopic spin current originates in the micro- 
scopic physics proposed in Ref. [g| : that superfluid elec- 
trons reside in circular orbits of radius 2A^ , their velocity 
given by Eq. (13), hence up and down spin electrons or- 
bit in opposite direction with orbital angular momentum 
antiparallel to their spin angular momentum and of the 
same magnitude, h/2. The state arises from the interac- 
tion of the electron with the positive ionic background, 
i.e. a spin-orbit interaction, as the wavefunction expands 
in the transition to the superconducting state to lower 
the kinetic energy of confinement of the antibonding elec- 
trons at the top of a nearly filled band. As shown in 
ref.JH], the expansion from a small orbit to a mesoscopic 
orbit of radius 2A l provides a dynamical explanation for 
the Meissner effect as well as for the dynamical origin of 
the spin current. 

Furthermore, the total mechanical angular momentum 
carried by each component of the spin current in the sur- 
face layer was found to be identical to the aggregate sum 
of the spin angular momenta of that spin component of 
the superfluid electrons in the entire volume. This fol- 
lows from the fact that the orbital angular momentum for 
each electron in its orbit of radius 2Xl is h/2. It implies 
that the entire superconductor is a single giant Cooper 
pair, with its giant spin current components reflecting 
as well as quenching the aggregate sum of the superfluid 
electron spins. 
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We also found that our equations predict that a sig- 
nificant amount of negative charge spills out of the su- 
perconductor. This effect was anticipated in our earlier 
work, and we believe it may play an important role in the 
proximity effect. It also clearly illustrates the tendency of 
superconductors to get rid of their excess negative charge 
(antibonding electrons at the top of the Fermi distribu- 
tion) predicted by the theory of hole superconductivity, 
which is also reflected in the predicted tunneling asym- 
metry of universal sign [29]] (larger tunneling current for 
negatively biased superconductor). 

The fact that in the proposed scenario spin-orbit cou- 
pling is an essential ingredient of superconductivity also 
provides a natural explanation for the mechanism by 
which the ionic lattice picks up the 'missing' mechani- 
cal angular momentum |26| when a metal is cooled into 
the superconducting state in the presence of a magnetic 
field. Conventional BCS-London theory neither provides 
an explanation for how the electrons in the Meissner cur- 
rent acquire their mechanical angular momentum, nor for 
how the ions acquire a compensating angular momentum 
in the opposite direction [26|. In our earlier attempt to ex- 
plain the angular momentum conservation puzzle [2q] we 
could not account for a complete Meissner effect precisely 
because the spin-orbit interaction mechanism to transfer 
angular momentum to the lattice was not included. For 
a 99% Meissner effect, we had to assume a charge expul- 
sion at least 3 orders of magnitude larger than the values 
discussed here (Eq. (27) of ref.JHj]). 

Finally, we discussed a new 'geometric' interpretation 
of type I and type II superconducting regimes and of the 
effect of disorder in superconductors based on the present 
model, and we proposed that simply connected type II 
superconductors cannot exist according to this theory. 

Note however that we have assumed at various points 
in this paper that the size of the superconductor is much 
larger than Al (e.g. Eqs. (2) and (3)). For superconduct- 
ing samples of dimensions comparable or smaller than 
the penetration depth we expect the fundamental equa- 
tions of Sect. VII to remain valid, but several features 
will change: in particular, the spin current at the surface 
will be smaller than the universal form Eq. (13), and 
the maximum electric field E m will be smaller than the 
value given by Eq. (33) Note also that our treatment 
assumed local electrodynamics and non-local corrections 
will be important in strongly type I superconductors. 

The hypothesis that electrons in superconductors re- 
side in orbits of diameter hundreds of lattice spacings 
was originally proposed by Slater 24j . We have shown 
that the radius of these 'Slater orbits' has to be pre- 
cisely 2Al to give rise to perfect diamagnetism. Indeed, 



as argued by Slater, large orbits provide a simple and 
compelling explanation for why superconductors cannot 
tolerate the presence of a magnetic field in their interior: 
it simply costs too much energy, proportional to the area 
subtended by these mesoscopic orbits. Type I supercon- 
ductors have only two ways to deal with this: either expel 
the magnetic field from their interior, paying the associ- 
ated electromagnetic energy price, or if the price is too 
high, become normal, collapsing the mesoscopic orbits to 
microscopic Landau orbits of radius kp 1 ~ lattice spac- 
ing and pay the associated kinetic energy price. Type 
II superconductors have a third way: enclose within the 
mesoscopic electron orbits tubes of magnetic flux where 
the system is normal. Type II superconductors are able 
to do this because the two members of a Cooper pairs 
have the centers of their 2Al orbits sufficiently close to 
each other that they can enclose within them the same 
flux tube. 

Experiments should be able to detect the existence of 
the predicted spin current near the surface of supercon- 
ductors. We have given numerical estimates for the mag- 
nitude of the spin current and other associated quantities 
for 'conventional' and other superconductors. The spin 
current density can be as large as 2 x 10 8 Amps /cm 2 , and 
should be detectable in inelastic neutron scattering ex- 
periments with very cold neutrons 30] . The excess charge 
near the surface as large as 1 electron per 10 6 ions and 
the internal electric field as large as 300, 000L/ cm should 
be experimentally detectable. It should be pointed out 
that the calculation of electric fields in the neighborhood 
of non-spherical samples discussed earlier [3l| needs to be 
modified in light of the results in this paper; this will be 
discussed elsewhere. Experimental verification of these 
predictions will support the basic tenets of the theory 
of hole superconductivity: that superconductivity origi- 
nates in the fundamental charge asymmetry of condensed 
matter, that it occurs when a metal has "too many elec- 
trons" (an almost full band, small de Broglie wavelength 
for the carriers at the Fermi energy, and with higher T c 
for anions [12]), and that it is driven by kinetic energy low- 
ering and involves dressed holes (antibonding electrons) 
turning into undressed electrons. Further development of 
the microscopic theory^, HH will be discussed in forth- 
coming work. 
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